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Abstract 

We introduce a model of proportional growth to explain the distribution P g {g) of business firm 
growth rates. The model predicts that P g (g) is exponential in the central part and depicts an 
asymptotic power-law behavior in the tails with an exponent £ = 3. Because of data limitations, 
previous studies in this field have been focusing exclusively on the Laplace shape of the body of 
the distribution. In this article, we test the model at different levels of aggregation in the economy, 
from products to firms to countries, and we find that the model's predictions agree with empirical 
growth distributions and size- variance relationships. 
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I. INTRODUCTION 

Gibrat , building upon the work of the astronomers Kapteyn |3j, assumed the ex- 
pected value of the growth rate of a business firm's size to be proportional to the current 
size of the firm, which is called "Law of Proportionate Effect" jj, 0|. Several models of 
proportional growth have been subsequently introduced in economics in order to explain the 
growth of business firms QDD Simon and co-authors extended Gibrat's 

model by introducing an entry process according to which the number of firms rise over time. 
In Simon's framework, the market consists of a sequence of many independent "opportuni- 
ties" which arise over time, each of size unity. Models in this tradition have been challenged 



by many researchers 



13, 



14. 
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16 



17| who found that the firm growth distribution is not 
Gaussian but displays a tent shape. 

Here we introduce a general framework that provides an unifying explanation for the 
growth of business firms based on the number and size distribution of their elementary 
constituent components 0, 0, 0, 0, • Specifically we present a model of 

proportional growth in both the number of units and their size and we draw some general 
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implications on the mechanisms which sustain business firm growth [?]. 
According to the model, the probability density function (PDF) of growth rates is Laplace 
in the center 1^] with power law tails 0,0] decaying as P g {g) ~ where ( = 3. 

Because of data limitations, previous studies in this field focus on the Laplace shape of 
the body of the distribution j3J|. Using a database on the size and growth of firms and 
products, we characterize the shape of the whole growth rate distribution. 

We test our model by analyzing different levels of aggregation of economic systems, 
from the "micro" level of products to the "macro" level of industrial sectors and national 
economies. We find that the model accurately predicts the shape of the PDF of growth rate 
at all levels of aggregation studied. 

II. THE THEORETICAL FRAMEWORK 



We model business firms as classes consisting of a random number of units. Accord- 
ing to this view, a firm is represented as the aggregation of its constituent units such 
as divisions businesses 0], or products 0]. Accordingly, on a different level of 
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coarse-graining, a class can represent a national economy composed by economic units such 
as firms. In this article we study the logarithm of the one-year growth rate of classes 
g = \og(S(t + 1)/S(t)) where S(t) and S(t + 1) are the sizes of classes in the year t and t + 1 
measured in monetary values (GDP for countries, sales for firms and products). Our model 
is illustrated in Fig. ^ Two key sets of assumptions in the model are (A) the number of 
units in a class grows in proportion to the existing number of units and (B) the size of each 
unit fluctuates in proportion to its size. 
The first set of assumptions is: 

(Al) Each class a consists of K a (t) number of units. At time t — 0, time step measured 
by year generally, there are N(0) classes consisting of n(0) total number of units. The 
initial average number of units in a class is thus n(0)/iV(0). 

(A2) At each time step a new unit is created. Thus the number of units at time t is 
n {t) = n(0)+t. 

(A3) With birth probability b, this new unit is assigned to a new class, so that the average 
number of classes at time t is N(t) = N(0) + bt. 

(A4) With probability 1 — b, a new unit is assigned to an existing class a with probability 
P a = (1 - b)K a (t)/n(t), so K a (t + 1) = K a (t) + 1. 

For simplicity, we do not consider the decrease of the number of units in a class. In 
reality, elementary units enter and exit. Because we are considering the case of a growing 
economy, it is legitimate to assume the entry rate being higher than the exit rate. On the 
average, the net entry rate of units can be simplified as a positive constant. In the model, 
the net entry rate of units is fixed at 1. Thus, at large t, it gives results equivalent to the 
ones that would have been obtained considering a value for the exit rate of units. 

Our goal is to find P(K), the probability distribution of the number of units in the classes 
at large t. This model in two limiting cases (i) 6 = 0, K a = 1 (a = 1, 2 . . . iV(0)) and (ii) 
6^0, N(Q) = 1, n(0) = 1 has exact analytical solutions P(K) = N(QVt(t/(t + N(0))) K (1 + 
0(1/*)) H Q and lim P(K) = (1 + b)T(K)T(2 + b)/V{K + 2 + b) Q respectively. 

t— »oo 

In the general case, the exact analytical solution is not known and we obtain a numerical 
solution by computer simulations and compare it with the approximate mean field solution. 

n 

(see, e.g., Chapter 6 of [35[ and Appendix A) 
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Our results are consistent with the exactly solvable limiting cases as well as with the 
empirical data on the number of products in the pharmaceutical firms and can summarized 
as follows. In the limit of large t, the distribution of K in the old classes that existed at 
i = 0_toane X ponen t ia,di Str ibutionQ 

PMK) = A A ' /f(f) 1 _ 1 « -^exp(-A-/A'( f )), (1) 

where A = 1 — 1/K(t) and K(t) is the average number of units in the old classes at time t, 
K(t) = [{n{0) + t)/n{0)] 1 ~ b -n{0) b /N(0). The distribution of units in the new classes created 
at t > converges to a power law with an exponential cutoff 

P new (K)~K-*f(K), (2) 

where (p ~ 2 + b for small b, and f(K) decays for K —>■ oo faster than P Q i^(K). The 
distribution of units in all classes is given by 

P{K) = Nwki p ^ + WWTTt < 3) 

The mean field approximation for P new (K) is given by 

PneUK) - ^±1^-1/(1^-1) f\- y yjh % (4) 
1 — ] K' 



where K' = K[n(0)/(n{0) + t)] 1 - 1 . 
The second set of assumptions is: 

(Bl) At time t, each class a has K a (t) units of size i = 1,2, ...K a (t) where i^Q, and 

£i > are independent random variables taken from the distributions P(K a ) and 
P^(Ci) respectively. P(K a ) is defined by Eq. (jBJ) and -P§(£i) is a given distribution with 
finite mean and standard deviation and ln£j has finite mean /i£ = (ln£j) and variance 
= ((ln£j) 2 ) — The size of a class is defined as S a (t) = J2i=i £i(X)- 

(B2) At time t + 1, the size of each unit is decreased or increased by a random factor 
rjiit) > so that 

&(* + !)= (5) 

where rji(t), the growth rate of unit i, is independent random variable taken from a 
distribution P v (r]i), which has a finite mean. We also assume that mr^ has finite mean 
fi v = (hirji) and variance V v = ((In r/j) 2 ) — /i 2 . 



The growth rate of each class is defined as 

g a ee log } ) = *og|>(t + 1) ~ log $>(*)■ (6) 

Here we neglect the influx of the new units, so K a = K a (t + 1) = K a (t). The resulting 
distribution of the growth rates of all classes is determined by 

oo 

P g (g) = J2P(K)P 9 (9\K), (7) 

where P{K) is the distribution of the number of units in the classes, computed in the 
previous stage of the model and P g (g\K) is the conditional distribution of growth rates of 
classes with given number of units determined by the distribution Pf(£) and P n {f])- 

The analytical solution of this model can be obtained only for certain limiting cases but 
a numerical solution can be easily computed for any set of assumptions. We investigate the 
model numerically and analytically (see Appendix B) and find: 

(1) The conditional distribution of the logarithmic growth rates P g (g\K) for the firms 
consisting of a fixed number of units converges to a Gaussian distribution for K — > oo: 



P g {g\K) » -p= exp {(g - g) 2 K/2V 9 ) , (8) 

where V g is a function of parameters of the distribution P^(£) and P v (t]), and g is mean 
logarithmic growth rate of a unit, g = (lnr/j). 

Thus the width of this distribution decreases as 1/ \fK. This result is consistent with 
the observation that large firms with many production units fluctuate less than small 



firms 



18, 



22, 



(2) For g 3> V v , the distribution P g (g) coincides with the distribution of the logarithms of 
the growth rates of the units: 

P g (g)^P v (]nr]). (9) 

In the case of power law distribution P{K) ~ which dramatically increases for 
K — > 1, the distribution P g (g) is dominated by the growth rates of classes consisting 
of a single unit K — 1, thus the distribution P g (g) practically coincides with P^ln^) 
for all g. Indeed, our empirical observations confirm this result. 



(3) If the distribution P(K) ~ K~ v , ip > 2 for K — > oo, as happens in the presence of the 
influx of new units b ^ 0, P s (<?) = C\ — C2\g\ 2ip ~ 3 , for g — > which in the limiting case 
6 — > 0, y? — > 2 gives the cusp P g (g) ~ Ci — C 2 |5 , | (Ci and C 2 are positive constants), 
similar to the behavior of the Laplace distribution Ph{g) ~ exp(— |p|C 2 ) for $ — > 0. 

(4) If the distribution P(K) weakly depends on K for K — > 1, the distribution of Pg(g') 
can be approximated by a power law of g: P g (g) ~ |g|~ 3 in wide range yVg~/ ' K{t) <C 
(7 <C \A^, where -ft'(i) is the average number of units in a class. This case is realized 
for b = 0, t — > oo when the distribution of P(K) is dominated by the exponential 
distribution and i^(t) — > oo as defined by Eq. (JTJ). In this particular case, P g (g) for 
g <C y/Vg can be approximated by 

(5) In the case in which the distribution P(K) is not dominated by one-unit classes but 
for K — > oo behaves as a power law, which is the result of the mean field solution for 
our model when t — ► oo, the resulting distribution -P 9 (g) has three regimes, P g (g) ~ 
Ci — C^lgl 2 ' 1 ' 5-3 for small g, -P 9 (g) ~ | C7 1 3 for intermediate g, and P g (g) ~ P(ln^) for 
<7 — > oo. The approximate solution of P g (g) is obtained by using Eq. (JSJ) for P g (g>|i^) 
for finite K, mean field solution Eq. (J3J) in the limit t — > oo for P(K) and replacing 
summation by integration in Eq. (JTJ): 

j | /*oo /*oo 

Pgi9) = T^b\y^vl Q M-y)y^ b dy J exp(-g 2 K/2V)K^-Th) dK . (ii) 

For 6 7^ the integral above can be expressed in elementary functions. In the 6^0 
case, Eq. (JTTJ) yields the main result 



2V^ 

~ ; ; (6 0) (12) 

9W VF+2^(|g| + v^T2^) 2 ' 



which combines the Laplace cusp for g — > and the power law decay | C7 1 3 for g — > oo. 
Note that due to replacement of summation by integration in Eq. (JTJ), the approxima- 
tion Eq. (JT2J) holds only for g < \fV r] . 

In Fig. we compare the distributions given by Eq. (|10j). the mean field approximation 
Eq. (jllj) for b = 0.1 and Eq. (JP2"j) for b 0. We find that all three distributions have very 
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similar tent shape behavior in the central part. In Fig. ^Bp we also compare the distribution 
Eq. (JT2j) with its asymptotic behaviors for g — > (Laplace cusp) and g — ► oo (power law), 
and find the crossover region between these two regimes. 

III. THE EMPIRICAL EVIDENCE 

To test our model, we analyze different levels of aggregation of economic systems, from 
the micro level of products to the macro level of industrial sectors and national economies. 

First, we analyze a new and unique database, the pharmaceutical industry database 
(PHID), which records sales figures of the 189,303 products commercialized by 7,184 phar- 
maceutical firms in 21 countries from 1994 to 2004, covering the whole size distribution for 
products and firms and monitoring the flows of entry and exit at both levels kindly pro- 
vided by the EPRIS program. Then, we study the growth rates of all U.S. publicly-traded 
firms from 1973 to 2004 in all industries, based on Security Exchange Commission filings 
(Compustat). Finally, at the macro level, we study the growth rates of the gross domestic 
product (GDP) of 195 countries from 1960 to 2004 (World Bank). 

Fig. El shows that the growth distributions of countries, firms, and products are well fitted 
by the distribution in Eq. ()12|) with different values of V g . Indeed, growth distributions at 
any level of aggregation depict marked departures from a Gaussian shape. Moreover, even 
if the P g (g) of GDP can be approximated by a Laplace distribution, the P g (g) of firms 
and products are clearly more leptokurtic than Laplace. Based on our model, the growth 
distribution is Laplace in the body, with power-law tails. In fact, Fig.0]show that the central 
body part of the growth rate distributions at any level of aggregation is well approximated 
by a double exponential fit. Fig. reveals that the asymptotic behaviors of g at any level 
of aggregation can be well fitted by power-law with an exponent ( = 3. 

Our analysis in Sec. II predicts that the power law regime of P g (g) may vary depending 
on the behavior of P(K) for K — > 1, and the distribution of the growth rates of units. In 
case of PHID, for which P(l) 3> P(2) 3> P(3) . . . the growth rate distribution of firms must 
be almost the same as the growth rate distribution of products, as we stated in Sec. II. Hence 
the power law wings of P g (g) for firms originate on the level of products. Because PHID 
does not contain information on the subunits of products we can not test our prediction 
directly, but we can hypothesize that the distribution of the product subunits (number 
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of customers or shipping ways) is less dominated by small K, but has a sufficiently wide 
power law regime due to the influx of new products. These rather plausible assumptions are 
sufficient to explain the shape of the distribution of the product growth rates, which is well 
described by Eq. 

The PHID database allows us to test the empirical conditional distribution P g (g\K) and 
the dependence of its variance a 2 on K, where K is the number of products. We find that 



u ~ 



K , which is significantly smaller than 1/yK behavior. This result does not imply 
correlations among product growth rates on the firm level 1 2 1 ] , but can be explained by 
the fact that for skewed distributions of product sizes -P§(£) characterized by large the 
convergence of P g (g\K) to its Gaussian limit Eq. (jSJ) is slow and the growth rates of the 
firms are determined by the growth of the few large products. Using the empirical values 
for the PHID fi^ = 3.44, = 5.13, ^ = 0.016, V v = 0.36 and assuming lognormality of the 
distributions P^) and P n {ff) we find that the behavior of a can be well approximated by a 
power law a ~ j{- - 20 for K < 10 3 . For this set of parameters, the convergence of P g (g\K) 
to a Gaussian distribution takes place only for K > 10 5 . This result is consistent with 
the observations of the power law relationship between firm size and growth rate variance 



reported earlier |13j, |18|, |19|, |3 



oi tne per 
|l3, 18, m 



IV. DISCUSSION 



Business firms grow by increasing their scale and scope. The scope of a firm is given 
by the number of its products. The scale of a firm is given by the size of its products. A 
firm like Microsoft gets few big products while Amazon sells a huge variety of goods, each 
of small size in terms of sales. In this article we argue that both mechanisms of growth are 
proportional. The number of products a firm can successfully launch is proportional to the 
number of products it has already commercialized. Once a product has been launched its 
success depends on the number of customers who buy it and the price they are willing to pay. 
To a large extent, if products are different enough, the success of a product is independent 
from other products commercialized by the same company. Hence, the sales of products 
can be modeled as independent stochastic processes. Moreover, sometimes, new products 
are commercialized by new companies. As a result, small companies with few products can 
experience sudden jerks of growth due to the successful launch of a new product. 
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In this article, we find that the empirical distribution of firm growth rates exhibits a 
central part which is distributed according to a Laplace distribution and power-law wings 
P g (g) ~ g^^ where ( = 3. If the distribution of number of units K is dominated by single 
unit classes, the tails of firm growth rate are primarily due to smaller firms composed of one 
or few products. The Laplace center of the distribution is shaped by big multiproduct firms. 
We find that the shape of the distribution of firm growth is almost the same in presence of a 
small entry rate and with zero entry. We also find that the model's predictions are accurate 
also in the case of product growth rates, which implies that products can be considered 
as composed by elementary sale units, which evolve according to a random multiplicative 
process Although there are several plausible explanations for the Laplace body of the 
distribution [3, LuL the power law decay of the tails has not previously been observed. 
We introduce a simple and general model that accounts for both the central part and the 
tails of the distribution. The shape of the business growth rate distribution is due to the 
proportional growth of both the number and the size of the constituent units in the class. 
This result holds in the case of an open economy (with entry of new firms) as well as in the 
case of a closed economy (with no entry of new firms). 
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Appendix A: The distribution of units in old and new classes 

Assume that at the beginning there are N(0) classes with n(0) units. Because at every 
time step one unit is added to the system and a new class is added with probability b, 
at moment t there are n(t) = n(0) + t units and N(t) = N(0) + bt classes, among which 
there are bt new classes with n new units and N(0) old classes with n u units, such that 
rioid + n new = n(0) + t. 

Because of the preferential attachment assumption, we have 

dTl new j , , 17"new / a 1 \ 
: = 0+(l— 0— r~. , All 

dt v ; n(0)+t' v ; 

dt V n(0)+t v ; 

Solving the second differential equation and taking into account initial condition n o ^(0) = 
n(0), we obtain 

n old (t) = (n(0)+t) 1 - b n{0) b . (A3) 

Analogously, the number of units at time t in the classes existing at time to is 

n e (M) = (ji(0)+t) 1 -\n{0)+t o ) b , 

where the subscript 'e' means "existing" . The average number of units in old classes is 

K(t) = = m+2± „ (0 )». (A4) 

v ; N(0) N(0) K J K ' 



It is known that 



for t —>■ oo the preferential attachment model converges to the 



exponential distribution: 

P app (K) w exp(-K/K(t))/K(t). (A5) 

Thus, we obtain 

Poid(K) w 7 7 . ^ . . „ exp f- / - 7 , (A6) 

V ; (n(0) +t) 1 - 6 n(0) b P V (^(0) +t) 1 - 6 n(0)V V ; 

and the part of P(K) of old classes is 

P oW (^)«P oW W^L. (A7) 

The number of units in the classes that appear at to is bdt and the number of these 
classes is bdt. Because the probability that a class captures a new unit is proportional to 

12 



the number of units it has already gotten at time t, the number of units in the classes that 
appear at time to is 

bdt 



(t ,t) = n e (t ,t) 



n(0)+t ' 

The average number of units in these classes is K(to,t) = n new (to,t)/bdt = (n(0) + 
t) 1 ~ b /(n(0) + to) 1 b - Assuming that the distribution of units in these classes is given by a 
continuous approximation in Eq. (|A5ft : 

Pne W {K) « — i_ exp (-K/K(t , t)) (A8) 
K{t ,t) 

Thus, their contribution to the total distribution is 

exp(-K/K(t ,t)) 



bdt 1 



N(0) + bt K(to,t) 
The contribution of all new classes to the distribution P(K) is 

p ~w K Nwrvt .[ xh) exp ( - K/K[t " t]) dt °- (A9) 

If we let y = K/K(t , t), then 

PnUK) ~ -±- K(-rh-i) ^1±L f K e -y y & dy . ( A 10) 

Note that Eq. (jA7|) and Eq. (|A10|) are not exact solutions but continuous approximations 
which assume K is a real number. Now we investigate the distribution in Eq. (j A 10ft . 

1. At fixed K when t — > oo, the low limit of integration in Eq. (|A10ft goes to zero and 
we have 



K 



Pne W {K) = i _ b I e- y y—dy. (All) 



As K — > oo, 



p^{k) = k-^(^) r(i + T l^). 



(A12) 



As K -> 0, 



P(K) = — — K^—^) _ = (A13) 



1-6 
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2. At fixed t when K — > oo, we use the partial integration to evaluate the incomplete T 
function: 



oo 

e~y y a dy = -e~ y y a \™ + a I e~ y y a ~ l dy » e~ x x 



Therefore, from Eq. ()A10|) we obtain 



new[ j iv(o) + &t i - 6 l ^ e y v " 



> n(0) + t J 



17 6X P I > ( Al4 ) 



iV(0) + bt 1-b K 1 \ \n(0)+t 
which always decays faster than Eq. (jA6|) . 

Appendix B: Calculation of the growth distribution of classes P(g) 

Let us assume both the size and growth of units (£j and 77, respectively) are distributed 
lognormally 

pfo) = -/=L= 1 exp (-(ln& - m f ) 2 /2^) , (A15) 

pfa) = —L= - exp (-(ln^ - m v ) 2 /2V v ) . (A16) 

If units grow according to a multiplicative process, the size of units = ^rji is distributed 
lognormally with = + V v and = + m v . 

The n th moment of the variable x distributed lognormally is given by 
f°° 1 x n 

fM x (n) = / —j= — dx exp (— (lnx — m) 2 /2V) = exp (nm x + n 2 V x /2) . (A17) 

Thus, its mean is \i x = = 1) = exp(m x . + V^/2) and its variance is a 2 = jj,2 — ^ 2 = 

ft 2 (exp(K)-l). 

Let us now find the distribution of g growth rate of classes. It is defined as 

9 = In = In f; £ - In f> (A18) 

^ {l> i=i i=i 

Here we neglect the influx of new units. According to the central limit theorem, the sum of 

K independent random variables with mean /i^ = and finite variance cx| is 

K 

J2ii = K^ + VKu K , (A19) 
i=i 
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where vk is the random variable with the distribution converging to Gaussian 

lim P(u K ) -»■ ; 1 exp i-vlllaf) . (A20) 
K -*°° Jto^* 

Because In fj, v = m v + V v /2 and In ^ — In fj,£ + In //,, we have 

£ = In 5(t + 1) - In S(t) = \n(Kn e ) + — ^ ln(A> f ) " K 



= ^ + ^Z!2*£. (A21) 

For large .fT the last term in Eq. (jA21j) is the difference of two Gaussian variables and that 
is a Gaussian variable itself. 

In fi£i = m^i + = In /zg + In /i^, (A22) 

where ln/i^ = + V v /2 is the average growth rate. To find the distribution of g we must 
find its mean and variance. In order to do this, we rewrite 

\K 



V 



and 

v k E£i(&-^e) 



Thus 



g = ra r 



= m v + Y± + E ^|^ f — ■ (A23) 

Because — A^/V the avera g e of each term in the sum is fx? ~^c^v = ®- The variance of 
each term in the sum is ((& rji) 2 ) — (2^f r)i fi v ) + ((,f ^1) where ^rji, £frji and are all lognormal 
independent random variables. Particularly, (^?7i) 2 is lognormal with V = 4V V + 4V^ and 
m = 2m v + 2m^ £?r)i is lognormal with V = 4V^ + V v and m = 2m^ + m^; is lognormal 
with V = AVj: and m = 2m ? . Using Eq. (lATYl) and Eq. (1X231) 

((&*7i) 2 ) = exp(2m„ + 2m ? + 2V„ + 2V ( ), (A24a) 
(C-?7i) = exp(m„ + 2m € + 2V C + V v /2), (A24b) 
(C 2 )=exp(2m c + 2\/ c ). (A24c) 
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Collecting all terms in Eqs. (jA24aHA24c|) together and using Eq. (|A23|) we can find the 
variance of g: 

2 K exp(2m ? + 2V 5 + 2m v + V v )(exp(V v ) - 1) 
° 9 ~ K 2 exp(2m € + V ( + 2m v + V v ) ' 

= iexp(V { )(exp(^)-l). (A25) 

Therefore, for large K, g has a Gaussian distribution 

/ > s v 7 ^ / (g-m) 2 K\ lk . 

where m = m,, + V^/2, V = exp(Vg)(exp(V^) — 1) and fj, v = exp(m v + V v /2). 

The distribution of the growth rate of the old classes can be found by Eq. (JJJ) in the 
text. In order to find a close form approximation, we replace the summation in Eq. (JZJ) 
by integration and replace the distributions P(K) by Eq. (jA6|) and P(g) by the Eq. (jA26|) 
assuming m = 0: 

1 f 00 1 K g 2 K 

2 vW V 2V y y ' V 7 

where is the average number of units in the old classes (see Eq. (|A4[) ). This distribution 
decays as 1/g 3 and thus does not have finite variance. In fact, we approximate the distri- 
bution of number of units in the old classes by a continuous function exp(—K/K(t))/K(t), 
while in reality it is a discrete distribution 



P old (K) = A ^— — 1 J , (A28) 

where A = exp(— 1/K(t)). The corrected distribution of growth rates is then given by the 
sum 

1 1 - A °° 

PoM « -T^rf^- E XK ^ exp(-g 2 K/2V). (A29) 



The slowest decaying term is 

1 



(1-A) exp{-g 2 /2V), (A30) 



V2nV 

which describes the behavior of the distribution when g — > oo. Thus there is a crossover 
from Eq. (IA"27|) to Eq. (lA"30l) when g « \/2K 
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For the new classes, when t — > oo the distribution of number of units is approximated by 



1 '* 



p ncw (K) « Y^i K ^~ b I y~ be ~ Vd y- ( A31 ) 







Again replacing summation in Eq. (J7J) in the text by integration and P(g\K) by Eq. (jA26j) 
and after the switching the order of integration we have: 

j | poo poo 

PneM « -— r, ^[W Jo eXp(_2/) ^ dy J exp( ~^ K/2V) # ( ~™ } dK - (A32) 

As g — > oo, we can evaluate the second integral in Eq. ()A32j) by partial integration: 

1 f°° 1 2V ill 

Pnew(g) ~ r / j== — y~T=b~* y 1 ^ exp(-y) exp(-yg 2 /2V) dy, 

1 ~ b Jo V2ttV g 2 

1 1 W 1 7-1. (A33) 



1 - b y^V 9 2 ^g 2 /2V + 1 £ 3 ' 

We compute the first derivative of the distribution ()A32j) by differentiating the integrand 
in the second integral with respect to g. The second integral converges as y — > 0, and 
we find the behavior of the derivative for g —>■ by the substitution K' = Kg 2 /(2V). As 
g — > 0, the derivative behaves as g • ^ 2 [ _ ( 3 / 2 ) +1 /( 1_b )] ~ g l /( l - b )- 2 ^ which means that the 
function itself behaves as C2 — Ci|<7| 2b /( 1 ~ b ) +1 , where C2 and G\ are positive constants. For 
small b this behavior is similar to the behavior of a Laplace distribution with variance V: 
exp(-V2\g\/VV)/V2V = 1/V2V - \g\/V. 

When b -> 0, Eq. (lA"32j) can be simplified: 

1 /"OO pK 

PneM\b->o « -7= / K- 3 / 2 exp(-K£ 2 /2VVK / expH/)j,dy, 



W V y/l + g y2V \g\/V2V+^g 2 /2V + l / 



Finally we find 



2V 

PncM)\^o « -.. , A34 



^2 

,3> 



which behaves for g — > as 1/v 2V — |p|/V and for (7 — > 00 as V/ (2<? ). Thus the distribution 
is well approximated by a Laplace distribution in the body with power-law tails. 

Because of the discrete nature of the distribution of the number of units, when g ^> \J2V 
the behavior for g — > 00 is dominated by const ■ exp(—g 2 /2V). 
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FIG. 1: Schematic representation of the model of proportional growth. At time t = 0, there 
are N(0) = 2 classes (□) and n(0) = 5 units (O) (Assumption Al). The area of each circle 
is proportional to the size £ of the unit, and the size of each class is the sum of the areas of 
its constituent units (see Assumption Bl). At the next time step, t = 1, a new unit is created 
(Assumption A2). With probability b the new unit is assigned to a new class (class 3 in this 
example) (Assumption A3). With probability 1 — 6 the new unit is assigned to an existing class 
with probability proportional to the number of units in the class (Assumption A4). In this example, 
a new unit is assigned to class 1 with probability 3/5 or to class 2 with probability 2/5. Finally, 
at each time step, each circle i grows or shrinks by a random factor rji (Assumption B2). 
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FIG. 2: (a) Comparison of three different approximations for the growth rate PDF, P g (g), given by 
Eq. (|10[). mean field approximation Eq. (jllj) for b = 0.1 and Eq. ()12l) . Each P g (g) shows similar tent 
shape behavior in the central part. We see there is little difference between the three cases, 6 = 
(no entry), b = 0.1 (with entry) and the mean field approximation. This means that entry of new 
classes (b > 0) does not perceptibly change the shape of P g (g). Note that we use K(t)/V g = 2.16 
for Eq. (|1()|) and V g = 1 for Eq. Q12j) . (b) The crossover of P g (g) given by Eq. (|12|) between the 
Laplace distribution in the center and power law in the tails. For small g, P g (g) follows a Laplace 
distribution P g {g) ~ exp(— \g\), and for large g, P g (g) asymptotically follows an inverse cubic power 
law P g (g) ~ g- 3 . 
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Scaled growth rate, (g-g)/V g 

FIG. 3: Empirical tests of Eq. (|12[) for the probability density function (PDF) P g {g) of growth rates 
rescaled by y/Vg. Shown are country GDP (0)> pharmaceutical firms (□), manufacturing firms 
(O), and pharmaceutical products (A). The shapes of P g (g) for all four levels of aggregation are 
well approximated by the PDF predicted by the model (dashed lines). Dashed lines are obtained 
based on Eq. (JI2J) with V g ps 4 x 1(T 4 for GDP, V g « 0.014 for pharmaceutical firms, V g ps 0.019 for 
manufacturing firms, and V g ~ 0.01 for products. After rescaling, the four PDFs can be fit by the 
same function. For clarity, the pharmaceutical firms are offset by a factor of 10 2 , manufacturing 
firms by a factor of 10 4 and the pharmaceutical products by a factor of 10 6 . Note that the data 
for pharmaceutical products extend from P g (g) = 1 to P g {g) ~ 10~ 4 and the mismatch in the tail 
parts is because P g {g) for large g is mainly determined by the logarithmic growth rates of units 
In?]. 
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Scaled growth rate, (g-g)/V g 

FIG. 4: Empirical tests of Eq. (|12|) for the central part in the PDF P(g) of growth rates rescaled 
by \/Vg- Shown are 4 symbols: country GDP (0)> pharmaceutical firms (□), manufacturing firms 
(O), and pharmaceutical products (A). The shape of central parts for all four levels of aggregation 
can be well fit by a Laplace distribution (dashed lines). Note that Laplace distribution can fit 
P g {g) only over a restricted range, from P g (g) = 1 to P g (g) ~ 10 _1 . 
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Scaled growth rate, (g-g)/V g 

FIG. 5: Empirical tests of Eq. (|12ft for the tail parts of the PDF of growth rates rescaled by y/Vg. 
The asymptotic behavior of g at any level of aggregation can be well approximated by power laws 
with exponents £ ~ 3 (dashed lines). The symbols are as follows: Country GDP (left tail: O) right 
tail: •), pharmaceutical firms (left tail: □, right tail: ■), manufacturing firms (left tail: O, right 
tail: ♦), pharmaceutical products (left tail: A, right tail: A). 
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